A statistical model with a standard Gamma distribution 
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We study a statistical model consisting of TV basic units which interact with each other by ex- 
changing a physical entity, according to a given microscopic random law, depending on a parameter 
A. We focus on the equilibrium or stationary distribution of the entity exchanged and verify through 
numerical fitting of the simulation data that the final form of the equilibrium distribution is that 
of a standard Gamma distribution. The model can be interpreted as a simple closed economy in 
which economic agents trade money and a saving criterion is fixed by the saving propensity A. Al- 
ternatively, from the nature of the equilibrium distribution, we show that the model can also be 
interpreted as a perfect gas at an effective temperature T(A), where particles exchange energy in a 
space with an effective dimension -D(A). 
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I. INTRODUCTION 

Statistical physicists like to understand how and why 
systems evolve from an initial toward an equilibrium 
macroscopic state. The equilibrium state, far from being 
just the "final state" of the dynamical evolution, actually 
reflects the details of the underlying dynamics. To write 
down the "microscopic equation" governing the dynamics 
of the evolution is a major goal. The various probabil- 
ity distributions, resulting from the different correspond- 
ing microscopic equations, have a relevant interest, in 
that they can be used to derive most of the macroscopic 
properties of the system. One of the foremost examples 
is the Maxwell-Boltzmann distribution for the velocities, 
which can be obtained as a solution of the equation which 
Boltzmann proposed for the evolution of the probability 
distribution for a dilute gas. 

One of the current challenges is writing down the "mi- 
croscopic equation" which would correspond to the cen- 
tury old Pareto law P| in Economics, which states that 
the higher end of the distribution of income f{x) follows 
a power-law 

where x is the income (money)and the exponent a has 
a value in the interval 1 to 2 [2, 0, 0, 01 • In the recent 
past, several studies have been made to investigate the 
characteristics of the real income distribution and provide 
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explanations (see Ref. |a] for a brief summary and more 
references) . 

It is our general aim to study a statistical model of 
closed economy which can be either solved exactly or sim- 
ulated numerically, analyze the microscopic equation and 
investigate what kind of macroscopic money distribution 
results |lllll|l3,|lll|l2. It would be of particular in- 
terest to see whether the study would give a clue as to 
why the Pareto law arises. 

In this paper, we study a statistical model consisting of 
N basic units which interact with each other by exchang- 
ing a physical entity x, according to a given microscopic 
law with one constant parameter A. We study the sta- 
tionary probability distributions /(x) for different values 
of the parameter A. We verify through numerical studies 
that the final form of the equilibrium distribution f{x) 
is that of a standard Gamma distribution. In Section 
II, we interpret the model as a simple closed economy in 
which economic agents trade money and a saving crite- 
rion is fixed by the saving propensity A. In Section III, 
using the nature of the equilibrium distribution, we show 
that the model can also be interpreted as a perfect gas 
at an effective temperature T(X), made up of particles 
exchanging energy in a space with an effective dimension 
r'(A). Finally in Section IV, we draw conclusions. 



II. THE MODEL ECONOMY 

We begin by considering a simple model of closed econ- 
omy, in which N agents can exchange money in pairs 
between themselves. All the agents can be initially as- 
signed the same money amount x, since this condition 
is not restrictive. Agents are then let to interact and, 
at every "time step", a pair (i,j) is randomly chosen 
and the transaction carried out. During the transaction, 
the agent money amounts Xi and Xj undergo a variation. 
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in which they are randomly reassigned between the two 
agents. The exchange law is such that the money is con- 
served during the transaction, 1 . G . X 2 j X J — X ^ I X ^ 

where x[ and x'^ are the money values after the transac- 
tion. This implies that at any time the average initial 
money x also represents the average money, (x) = x . 



(1 — A). The new law, which replaces Eqs. Ql, is 

a;- = Xxi + e{l - X){xi + Xj) , 

x'j — Xxj + (1 - - )^}{xi + Xj) . 

These equations can also be rewritten as follows. 



(3) 



Case with A = 



The exchanges are made according to the following law: 



e (X; 

a- 



+ Xj) , 

e){xi + Xj) 



(1) 



where e is a random number uniformly distributed in the 
interval (0,1). It can be noticed that, in this model, 
agents have no debts after the transaction, i. e. they are 
always left with a money amount a; > or, equivalently, 
that Xi is a positive definite quantity, if a; > 0. 

It can be shown that, as a consequence of the conserva- 
tion of money, the_system relaxes toward a Gibbs money 
distribution m \~ 



(2) 



where (x) represents the average money. This means 
that, after the relaxation, most of the agents have a very 
small amount of money, while the number of very rich 
agents is exponentially small. In other words, for a given 
x' > 0, the number of agents with x > x' , as well as the 
total amount of money they own, exponentially decreases 
with x' . The equilibrium state represented by the Gibbs 
distribution Q has been shown to be robust, in that it is 
reached independently of the initial conditions and also 
in models with multi-agent transactions. 

We have re-obtained the exact Gibbs solution for the 
case A = by numerical simulations of a system with 
A'' = 500 agents, each agent having initially a money 
amount a; = 1. The system was evolved for 10^ time 
steps - i. e. transactions - in order to reach equilibrium, 
and the final equilibrium distributions was averaged over 
10^ different runs. Figure ^ shows that the numerical 
results (open circles for A = 0) are in good agreement 
with the Gibbs distribution (continuous line). 



B. 



Case with A > 



We now introduce a saving criterion in the same model 
economy, by assuming that the saving propensity, which 
represents the fraction of money saved before carrying 
out the transaction, is non-zero, i. e. A > ItI,!^. Con- 
servation of money still holds, Xi + Xj = x'^ + x'p but the 
money which can be re-assigned in a transaction between 
the i-th and the j-th agent has now decreased by a factor 



X 



Xa 



■ = Xi+ /^x , 



"3 — ^] 

Ax = {l-\)[ex.j-{l~€)xi\ 



A 



X 



(4) 



in which money conservation is manifest. 

We studied the equilibrium distribution of this model 
through numerical simulations, for various values of A, for 
N — 500 agents, again each agent having money x = 1 in 
the initial state. In each simulation a sufficient number 
of transactions, as far as 10^, depending on the value 
of A, was used in order to reach equilibrium. The final 
equilibrium distributions, for a given A, were obtained 
by averaging over 10^ different runs. The numerical data 
are shown in Fig. (cases A ^ 0). 




FIG. 1: Equilibrium money distributions for different val- 
ues of the saving propensity A, in the closed economy model 
defined by Eqs. 10 and Q. The continuous curves are the 
fitting functions, defined in Eq. ||SJ with the values of n gov- 
erned by Eq. Q. Note that for the simulation (x) = 1. 

We found an analytic form for the equilibrium distri- 
bution, for a given A (0 < A < 1), which turns out to fit 
extremely well all data ■ The function is conveniently 
expressed in terms of the parameter 



n{\) 



1 



3A 

1-A 



(5) 



This particular form of n(A) was suggested by a mechan- 
ical analogy, discussed in Sec. IIIII between the closed 
economy model with TV agents and the dynamics of a gas 
of N interacting particles. Then the money distributions. 
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FIG. 2: Same quantities as in Fig. but on a linear- 
logarithmic scale. Note that for the simulation {x) — 1. 



for arbitrary values of A, are well fitted by the function 

fnix) = a„a;""^ exp (-nx/ (x)) , 

where n is defined in Eq. jsj and the prefactor a„, where 
r(n) is the Gamma function, is fixed by the normaliza- 
tion condition dxfn(x) — 1. 

The fitting curves for the distribution (continuous 
curves) are compared with the numerical data in Fig. ^ 
The fitting describes the distribution also at large values 
of shown by the linear-logarithmic plots in Fig.|21 In 
Fig. 121 the numerical values of the parameters n{X) and 
On (A) obtained directly from fitting the data (shown as 
dots) are compared with the respective fitting functions 
and (shown as continuous curves). 

By introducing the rescaled variable 

e = nx/{x) , (7) 
the probability distribution © can be rewritten as 
1 



fn{x) — T^/ N 

n 1 (n) 



r-'exp(-e)EE7„(C) , (8) 



where 7n(0 is the standard Gamma distribution 0,^3- 
The cumulative distribution for 7n(0 is the incomplete 
Gamma function r(^, n) = J^°° d^' 7„(^')j 



d r(e,n) 

r(n) 



(9) 



The Gibbs distribution (|2Jl is a special case for n — 1. 
The term {x) /n, on the LHS of Eq. is just the scahng 
factor appearing when the change of variable, from ^ to 
X, is made in the last equation above in order to obtain 
the distribution fn{x) for the variable x. 




FIG. 3: The parameters n (top) and a„ (bottom) versus A, 
obtained from numerical data (dots) and the corresponding 
analytical formulas (continuous curves) given by Eqs. Q and 
@, respectively. Note that for the simulation {x} = 1. 
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FIG. 4: Variation of the mode with the parameter A. Dots 
represent the numerical data, while the dashed curve was ob- 
tained theoretically. 



We notice that, with respect to the Gibbs distribu- 
tion (0), the distribution defined by Eqs. ^ contains the 
power x"^^ and the factor n in the exponential, which 
qualitatively change the distribution shape. First, they 
lead to a mode, Xrm different from zero. The mode 
is shown as a function of the parameter A in Fig. 01 
where the dashed curve represents the theoretical predic- 
tion that the mode Xm — 3A/ (1 -I- 2A) obtainable from 
Eq. ®. Secondly, the presence of the factor n is rele- 
vant for the mechanical analogy, considered in detail in 
section imi Finally, in the limit A — > 1 (i. e. n — > oo), the 
distribution fn{x) tends to a Dirac (5- function, peaked 
around the average value (x). A qualitative picture of 



4 



the evolution of the shapes of the /„(x)'s, for A going 
from zero to unity, is obtained by inspection of the var- 
ious curves in Fig. A more rigorous derivation of the 
asymptotic distribution for A — > 1 can be made by study- 
ing the characteristic function </)(g) The Gamma 
distribution 71 (^) for the dimensionless variable £, has a 
characteristic function <j)i{q) = /g °° exp(i(7^)7i (^) — 
1/(1 — zg). The characteristic function of the Gibbs dis- 
tribution Ji{x) in Eq. (0) is obtained by rescaling q by 
the constant factor xj^— (x), 



the reduced variable ^ = Dx/2 (x) 



D' 



'/'I (9) = (1 {x)y 



(10) 



The characteristic function of the generic Gamma dis- 
tribution TaXC) is simply given by the n-th power of 
Ml) HI 111, Mq) = 1/(1 -ig)". Analogously, the 
corresponding characteristic function of fn{x) Eq. ((HJ, is 
obtained by scaling q hy x/S, — (x) /n, 

0„(g) = (l-zg(x)/n)-" . (11) 

Thus, in the limit n — > cxd (A — > 1), one obtains 

0„(g) exp{iq{x)) . (12) 

The corresponding distribution is obtained by transform- 
ing back the characteristic function, i. e. 

/+00 
dq exp{~iqx)(l)n{q) — > S{x ~ {x)) . 
-00 

(13) 

This limit shows that a large saving criterion leads to 
a final state in which economic agents tend to have the 
same amount of money and, in the limit of A — > 1, they 
all have the same amount (x). 



III. THE GAS MODEL 

The equilibrium distributions ^ can also be inter- 
preted as the Gibbs distribution of the energy x, for a gas 
at temperature T = (x) /ks- This establishes a link be- 
tween the type of closed economy models considered here 
and statistical systems, suggesting a re-interpretation of 
the economy model in terms of a mechanical system of in- 
teracting particles. The introduction of a saving param- 
eter A > changes the shape of the Gibbs distribution 
into that of a Gamma distribution, but the correspon- 
dence with a mechanical system is lost only apparently. 
In fact, the Gibbs distribution (0) can represent the dis- 
tribution of kinetic energy x only in D — 2 dimensions, 
when its average value is given by {x)2 — 2{kBT /2). In 
all other cases (Z) 7^ 2), it is easy to show, starting from 
the Maxwell-Boltzmann distribution for the velocity in 
D dimensions, that the equilibrium kinetic energy distri- 
bution f{x) coincides, apart from a normalization factor, 
with the Gamma-distribution 7n(^) with n = D/2 for 



D 



D/2 



r(#) 



cxp 



Da 



2{x) 



D 



D{x), 



DkeT 



(14) 



where {x) j-, represents the average value of kinetic energy 
in D dimensions. The analogy between the factor D in 
the argument of the exponential function in Eq. H14|l and 
the analogous factor n in Eq. JHJl is to be noticed. The 
main difference is that, while D is an integer number by 
hypothesis, the parameter n(A) can assume in general 
any real values larger than or equal to one. 

In Eq. 1)14(1 temperature appears implicitly as T = 
2 {x) /ksD. This suggests that also in the closed econ- 
omy model considered above the effective temperature of 
the system should be defined as (x) /n, rather than (x). 
This is a natural consequence of the fact that the average 
value of kinetic energy in D dimensions is proportional 
to D, due to the equipartition theorem, and that an es- 
timate of the amplitude of thermal fluctuations, which is 
independent of its effective dimension, can be obtained 
from the ratio {x)j^ /D. 

Direct comparison between Eqs. (|14l) and © leads to 
a formal but exact analogy, between money in the closed 
economy model considered above, with N agents, saving 
propensity < A < 1, and given average money (x), 
on one hand, and kinetic energy in an ensemble of N 
particles in D dimensions at temperature T, on the other, 
if the effective dimension and temperature are defined as 



D{X) = 2n{X) = 



T(A) = 



n(A) 



2(1 + 2A) 
1- A ' 
1-A 
1 + 2A ' 



x) 



(15) 



respectively. This equivalence can be qualitatively under- 
stood in terms of the underlying microscopic dynamics by 
considering the example of a fluid of interacting particles. 
In one dimension, particles undergo head-on collisions, in 
which they can exchange the total amount of energy they 
have. In an arbitrary (large) number of dimensions, how- 
ever, this is not possible for purely kinematic reasons and 
only a fraction of the total energy is actually released or 
gained on average in a collision. Since the equipartition 
theorem implies that on average kinetic energy is equally 
shared among the D dimensions, one can expect that, 
during a collision, only a fraction ^ 1/D of the total 
energy is exchanged (and that a corresponding fraction 
A ~ 1 - l/D is "saved"). This estimate ~ l/Z? of the ex- 
changed energy is to be compared with the expression for 
the fraction of exchanged money obtained from Eq. Q 
using n — D/2, namely 1 — X ^ 3/ {D /2 + 2), which was 
in fact found starting the fitting of the numerical data 
from a function prototype of a form similar to 1/D. 
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IV. CONCLUSIONS 

We have studied a statistical model, which can be in- 
terpreted as a generalization of the simple closed econ- 
omy model, in which a random reassignment of the total 
agent money a;, involved in the transaction, takes place. 
The generalized model is characterized by N agents car- 
rying out transactions according to a saving criterion, 
determined quantitatively through a saving propensity 
A > 0. Alternatively, it can be considered as represent- 
ing a gas of N interacting particles which on average 
exchange only a fraction of their kinetic energy x, dur- 
ing a collision. We have shown the existence of such an 
analogy, by empirically obtaining the corresponding an- 
alytical solution fn{x) for the equilibrium distribution 
from numerical data. 

In both cases the equilibrium distribution can be writ- 
ten as a Gamma distribution 7n(^), where the reduced 
variable is given by ^ = nx/ (x)^ and {x)^ represents 
the the average value of x, given by (x}„ = n (x}^. The 
equivalence is represented by n = n(A) = 1-1- 3A/(1 — A) 
being a function of the saving propensity on one hand 
((x)j^ is the average value for A = 0), and by n = D/2 
being just the half of the number of dimensions, on the 
other hand ((x)^ being in this case the average value in 
two dimensions). 

The fact that we obtain basically the same equilib- 
rium distribution, characterizing the kinetic energy a gas 
of particles, suggests some general considerations about 



closed-economy models. The mechanical analogy illus- 
trated above can be addressed to the fact that the sys- 
tem is described statistically by a micro-canonical en- 
semble, just as a closed mechanical system, in which the 
exchanged quantity is conserved. Thus, a saving propen- 
sity larger than zero or any other change in the micro- 
scopic law can be expected to lead to a different shape 
of the equilibrium distribution, as shown in the present 
work, in which e.g. the effective number of dimensions 
and temperature may be different. However, the simple 
fact that the money is conserved implies that one cannot 
obtain an arbitrary distribution but, rather, only equilib- 
rium distributions related directly to the micro-canonical 
ensemble. It is suggested that open systems may be char- 
acterized by an equilibrium or stationary state character- 
ized by a type of power-law distribution. 
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